Raman Sideband Cooling in presence of Multiple Decay Channels 
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We have investigated the efBciency of pulsed Raman sideband cooling in the presence of multiple 
decay and excitation channels. By applying sum rules we identify parameter regimes in which 
multiple scattering of photons can be described by an effective wave vector. Using this method we 
determine the rate of heating caused by optical pumping inside and outside the Lamb-Dicke regime. 
On this basis we discuss also the efficiency of a recently proposed scheme for ground-state cooling 
outside the Lamb-Dicke regime [G. Morigi, J.f. Cirac, M. Lewenstein, and P. ZoUer, Europhys. Lett. 
39, 13 (1997)]. 
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, Laser-cooling allows to cool ions and atoms to very low temperatures. For this purpose, the full knowledge of 
■ the effects of the various physical parameters determining the cooling process is very important. Among the various 
schemes, Raman sideband cooling has been demonstrated to be a very successful technique for preparing atoms in 
the ground state of a harmonic potential . This cooling method exploits two stable or metastable atomic internal 
levels, which we call \g) and |e), connected by dipole transitions to a common excited state \r). The transitions are 
, usually driven by alternating pulses. A typical sequence alternates a coherent pulse, in which the atom is coherently 
transferred from \g) to |e) via a properly designed Raman pulse, with a re-pumping pulse, in which the atom is 
incoherently re-scattered to \g) by means of a laser resonant with |e) \r). A change of the motional state during the 
repumping is a process of higher order in the ratio ur/v of the recoil frequency wr = Hk'^ /2m and the trap frequency 
v, with m being the mass of the atom and k the wave vector of the one-photon transition. In the Lamb-Dicke regime, 
where dr ^ the probability for a change of the motional state is negligible and therefore, on the average, the 
' system is cooled at a rate of one phonon of energy hv per cooling cycle. Since there is a finite probability for the atom 
C^, to be returned to the state |e) instead of being repumped, a number of incoherent scattering events may be required 
before the atom is finally scattered into \g), which significantly increases the motional energy at the end of the optical 
^ [ pumping, reducing the cooling efficiency. Furthermore, since two and three level schemes are realized using Zeeman 
^ or hyperfine substates, decays from |r) into other electronic substates can occur, leading to additional heating. 

In this work we quantify the effect of a finite branching ratio in pulsed Raman sideband cooling by calculating 
• • . the average shift and diffusion of the vibrational energy distribution at the end of an incoherent pumping pulse. It 
should be pointed out that theoretical studies on laser-cooling for multilevel ions exist, which systematically include 
the branching ratio in their treatments l^-^ . Those studies have focussed on the Lamb-Dicke regime and on certain 
^ ' cooling schemes. Here, we single out the effect of the branching ratio on cooling for an arbitrary ratio ojr/v by applying 
5^ , sum rules. Hence, we infer the cooling efficiency in the Lamb-Dicke regime and we discuss the result outside the Lamb- 
Dicke regime in connection with the proposal in In particular, we show that in some parameter ranges the average 
effect of the multiple photon scattering can be described with an effective wave vector k^ff for the "effective" two-level 
transition \e) \g) 

This article is organized as follows. In Section 2 we introduce the model for the evolution of a trapped ion during 
the repumping pulse in a Raman transition, and we evaluate the average shift and variance of the ion energy at the 
end of the pulse. In Section 3 we extend our analysis to cases where the channels of decay are multiple. In Section 4 
we draw some conclusions, and in the Appendix we report the details of our calculations. 



II. MODEL 

We consider a three level atom as in Fig. |^, whose internal levels are a ground state |g), stable or metastable state 
|e) and excited state |r) of radiative width 7; |g) — > |r), |e) — > |r) are dipole transitions, with respective probabilities 
of decay pg, Pe, where Pg + Pe — 1. A laser resonantly drives the transition |e) — + |r) with Rabi frequency f2e. In the 
following we assume the wave vectors for both transitions to be equal to k, which is a good approximation if, e.g., |e) 
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and |g) are hyperfine components of the ground state. We study the ion motion in one-dimension. 
The master equation for the atomic density matrix is written as {h = 1): 

^P3 = -i [Ho + V, P3] + Lp3, (1) 

where Hq has the form: 

Ho = 6\e){e\ + iya+a. (2) 

Here, S is the detuning of the laser on the |e) — > \r) transition, which we take to be zero, and v is the frequency 
of the harmonic oscillator which traps the ion along the x-direction, with a, a+ annihilation and creation operator, 
respectively. The interaction of the ion with the laser light is described in the dipole approximation by the operator 
V: 

V^^{Ate^'^+h.c.), (3) 

with Aj' = (with j ~ e,g) dipole raising operator, Aj its adjoint, and x the position of the atom. In writing 

(0) have applied the Rotating Wave Approximation and we have moved to the inertial frame rotating at the 
laser frequency. Finally, the relaxation super-operator has the form 

iP3 = -J(|r)(r|p3+P3|r)(r|) 



+ f du7V(«)A7e-^^"^p3e^'="^A+, 



,2^ 



where N{u) is the dipole pattern of the spontaneous emission, which we take 7V(u) = 3/8(1 + u^). 
In the limit fig ^ 7 we can eliminate the excited state \r) in second order perturbation theory |Q, and reduce the 
three-level scheme to a two level one, with excited state |e) and linewidth 7e = In the limit v the 

master equation for the density matrix p, projection of p^ on the subspace {|e), \g)}^ can be rewritten as 



with iJcff effective Hamiltonian 



p^^i [iJeftP - pH+f^] + 7e [JeP + Jgp] , (4) 



ffcff = i^o-»y|e)(e|, (5) 



and with J^p, JgP jump operators, defined as: 

JjP = P] i'^je [J'p] <^ej) with j = g, e, (6) 

where cry = and where 

J'p = j dMiV(u)e-*'=(i+")^pe*'=(^+")^. (7) 
The solution of Eq. (^ can be written as follows ^ : 

p{t) = S{t)p{Q) + 7e / dtiS{t - ti)JS{ti)p{0) + ... (8) 

JQ 

+ le I dtl r dh... /" ' dtr,S{t-h)JS{ti-t2)Je-JeS{tn)p{(i) + ... . 







with J = Je + Jg, and S(t) is the propagator for the effective Hamiltonian: 

S{t)p{0) = e-*^=«*p(0)e'^"«'. (9) 



In Eq. (g) the successive contributions to the multiple scattering event are singled out: The first term on the RHS 
corresponds to the case in which at time t no spontaneous decay has occurred. The second term describes a single 



2 



scattering event, and the n-th term n — 1 scattering events. The trace of each term corresponds to the probabihty 
associated with each event, and we can thus interpret Eq. (||) as the sum over ah the possible paths of the scattering 
event weighted by their respective probabihties. At i — > cx), p{t) ps , the atom is in \g) and (e|ps|e) = 0. For a 
pulse of duration t ^ 1 j^f, we can replace i by oo in the integrals of Eq. (H) and assume that the atom has been 
scattered into \g) at the end of the pulse. Now, each term on the RHS of Eq. (||) corresponds to the path associated 
with a certain number of scattering events into |e) before the atom is finally scattered into \g). Through we can 
evaluate the shift and the variance of the energy distribution at the end of the repumping pulse, which are defined as: 

A^ = Tr{(7?„,ec-So)ps}, (10) 

OE = ^Jti-{{H,^,, - Eo- AEf ps}, (11) 

where i/mcc — va^a and Eq is the initial motional energy of the atom. 



A. Evaluation of the average shift and diffusion 

For simplifying the form of the discussion presented below, we rewrite the operator J' as follows: 

J' p — Jp+ Jp, 

where J, J are defined as: 

Jp = Ei\i){i\[Jp]\i){il 
Jp^EiEi,,i,^i\i){i\[Jp]\h){hl 



(12) 

(13) 
(14) 



and where {\l)} is the basis of eigenstates of the harmonic oscillator. For p{0) = |e)(e| (g) /i(0), with ^(0) initial 
distribution over the motional states, and according to Eq. (@) the steady state distribution has the form: 



^^00 = ^ PgPT-'ry{o) + F{j, j)/i(o). 



(15) 



where /Joo = (dlpsld) is the final distribution over the motional states. The first term in the RHS of ( p^ is the sum 
over all paths from |e) into \g), where after each jump the density operator is diagonal in the basis whereas the 

second term contains all other paths. These latter terms can be neglected |^], and for p{0) — \n){n\ the following 
relation holds: 



{S\Poo\s) 



£ J"|n)(r 



Dn{s). 



(16) 



Here, Dn{s) is the probability for the atom to be found in the state \g, s) at i — > oo, given the initial state |e,n) at 
t = 0. Using the explicit form (|l3|) of J in (|6|), -D„(s) has the form: 



m=0 fc=0 ^1 

oo oc 

■ J2 ■■ |(s|e''=''°(i+"i)(''*+'')|fci)p...|(fc™_i|e*'^'"°(i+"'")(''*+''V> 



(17) 



fei=0 fc„ 



where we have used the relation x = ao(a^ +a), with oq = \/\l2mv size of the ground state of the harmonic oscillator. 
Substituting ( p7| ) into Eqs. (p^, (|Tl|), and applying the commutation properties of 0,0+ [see the Appendix], we find: 



7 1 



2 2 



5 1 - Pe 



1 



5^"" ■ ''\-p, V 5 ; 49(l-pe)2 
where t] — kao = \/ ljJr/v is the Lamb-Dicke parameter. 



(18) 
(19) 
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B. Discussion 



Equation ( p8| ) represents the average shift to the vibrational energy at the end of the repumping pulse. For pe — 
it corresponds to the average recoil energy cj^ associated with one incoherent Raman scattering into \g). In this case, 
the second term in the RHS of Eq. (^9|) vanishes, and Eqs. (p^), ( p^ ) describe the scattering of one photon of wave 
vector k' — ^Jljhk on the effective two- level transition |e) \g). Similarly for > an effective wave vector /ceff 
can be defined for the incoherent scattering on the two-level transition |e) |g), which has the form 

fc' /7 fc 

VI -Pe V 5 VI -Pe 

Thus, fceff describes the average mechanical effect on the ion resulting from the multiple scattering of photons during 
the repumping pulse in a Raman transition with branching ratio (1 — Pe)/pe'- This description is valid in the limit in 
which we may neglect the second term in the RHS of ([l9|), i.e. for pe and/or 77 sufhciently small. In Fig. 2 the first 
term of RHS of Eq. (|l^) is compared with the complete expression for n — 0, for different values of the Lamb-Dicke 
parameter and as a function oipe. Here, we see that fcoff characterizes the scattering process for almost any branching 
ratio in the Lamb-Dicke regime, whereas for 77 = 0.6 an appreciable difference is already visible at pe = 0.2. 
From ( pO[ ) we can define the effective Lamb-Dicke parameter -q^s = fceff^o describing an incoherent scattering into 
the state \g). This parameter provides an immediate estimate of the effect of the branching ratio on cooling. For 
r|^/n ^ 1, if rjcsy/n ^ 1 the system is still in the Lamb-Dicke regime once it has been finally scattered into \g). 
Furthermore, the coarse-grained dynamics of the system can be described by a rate equation for the motional states 
\n) projected onto \g), where the rate of cooling (heating) is the real part of the sum of two terms: one corresponding 
to the component of the fluctuation spectrum of the dipolc force at frequency v {—v), the other to the diffusion 
coefficient due to spontaneous emission from the excited state [^[^. This latter term is proportional to the squared 
Lamb-Dicke parameter for the incoherent scattering, and thus in our case to yy^g. From the well-known solution of the 
rate equation [ p!o| , the diffusion term affects the steady state average vibrational number (n), which is proportional 
to the diffusion coefficient. 

Outside the Lamb-Dicke regime, when lor is comparable to, or larger than, i>, there are no estabilished ground-state 
laser-cooling techniques for trapped atoms. Here, we discuss our result in connection to the proposal in Q. There, 
a cooling scheme similar to Raman sideband cooling has been presented, where pulses which pump the atoms to the 
ground state alternate with pulses confining the atoms to a limited region of motional energy. These confinement 
pulses have two-photon detuning 5c to the red of the two-photon resonance frequency, where Sc ~ w^. Then, the 
presence of a branching ratio must be taken into account by choosing dc ~ AE. In this regime, pulses which efficiently 
counteract the average kick AE can be designed, provided that the following condition is fulfilled: 

fc-i^ > fccff, (21) 

where k^°^ is the projection on x of the two-photon wave vector of the coherent pulse. For two counterpropagating 
beams parallel to x, fc™*^ = 2\k\ and (|l|) is fulfilled for Pe < 13/20, i.e. up to branching ratios Pe/Pg ~ 2. Finally, 
outside the Lamb-Dicke regime the second term in the RHS of Eq. ( [l9| ) cannot be neglected. Hence, the diffusion is 
larger, and the efficiency of cooling may decrease dramatically as Pe increases. 



III. EXTENSION TO MULTI-LEVEL SCHEMES 



In the following, we show that the average heating associated with the repumping pulse in multilevel-schemes can 
be described in the same way as discussed in the previous sections. 

Let us consider the level-scheme of Fig. 3(a), where we have added to the scheme of Fig. 1 a further channel of decay 
from |r) into the stable or metastable state |1), with probability of decay pi such that Pi + Pe + Pg = 1, where p'^, 
p'g are the probability of decay onto |e), respectively. A laser resonantly drives the transition |1) ^ \r) with Rabi 
frequency For 7 the state |r) can be adiabatically eliminated from the equations of motion. In this 

limit the Master Equation aquires the form 

±p = -^[H,sp-pH+] (22) 

+P'eiJeP + Pll'JlP + p'gl'JgP, 

where 7' = 7e + 7i, with 7^ = f^j/7. The effective Hamiltonian is now: 
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i?efr-ffo-*f |e)(e|-z^|l)(l|, (23) 

and the jump operators have the form: 

■hp = XI ("^ ^ ^'^^^ 
with i = 1, e, g. The solution at t — > oo can be written as: 

oo 

Mcc^- Ep;(p^+Pi)"-^JXO). (25) 

m— 1 

Hence, the shift and variance have the form evaluated in Eqs. (p^), ( pj[ ) where now the probability Pg, pe are defined 
as Pe = p'e + Pi, Pg = p'g {pg +J3e = 1)- In a similar way we have evaluated these quantities for schemes like the 
one shown in fig. 3(b), where a second excited state |2) is coupled to |e) via the same recycling laser tuned on the 
transition |1) \r). For simplifying the treatment, we assume that a fourth laser resonantly drives the transition 
|1) — *■ |2) with Rabi frequency Vl (grey arrow in fig. 3(b)). Thus, for low saturation Eq. ( ^ ) describes the dynamics, 
where now 7^ — + 7,^'' {i = e, 1), with 7P'' being the rate of scattering through the excited state \i) [j = r,2). 

Assuming that Q, is such that ^If'^ = 7i'^V7i^'' = '^1 the solution in Eqs. (|l8|), (19) applies to this case too, where 
now Pf, is defined as: 

Pe = (p'e+Pl)T^ + -^, (26) 

1 + a 1 + a 

and the probability Pg of decaying into \g) is pg — 1 — p,,. 

The result ( |2^ ) shows that the total heating is minimum for a ^ 1, which can be obtained by choosing properly the 
laser intensity of the repumping lasers, or simply by removing degeneracies in the Zeeman multiplet, for example with 
the help of a magnetic field. 



IV. CONCLUSIONS 



We have studied the motional heating associated with a finite branching ratio and in the presence of multiple decay 
and excitation channels at the end of a repumping pulse in Raman sideband cooling. The first and second moments 
of the final energy distribution has been evaluated analytically, and the effect of the branching ratio has been singled 
out. We have shown that in a certain range of parameters the diffusion can be described with an effective wave vector 
fcoff, corresponding to an effective Lamb-Dicke parameter rjes for the incoherent scattering on the two- level transition 
|e) — > \g). Finally, on the basis of this result we have discussed the efficiency of Raman sideband cooling and of a 
recent proposal of ground-state cooling outside the Lamb-Dicke regime 

Analogous sum rules and considerations can be applied to Raman cooling for free atoms In that case the 

calculations are much simpler, since the total momentum of radiation and atom is a conserved quantity in the 
scattering event. 

In general, these results can be applied to cooling schemes in multilevel atoms. 
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VI. APPENDIX 



Using (p^, we rewrite ( p^ ) and (11) as: 
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oo 

4 = E 



n-AE/i^f{s\ J"»(n| |s), 
J"|n)(n| 



s . 



where we have introduced the quantities 

oo 

oo 
s=0 

Using Eq. (0), Eq. (|^) is rewritten as: 

oo „i „i 

^™ = rf"' E - ") / ^"1- / du„,iV(ui)...iV(u„) 

OO OO 

• E ••■ E l(^i|e"'^^+"'^'"^+''^|fe>P |(fc„|e"'(i+"'")('^*+"V)P- 



(27) 
(28) 

(29) 
(30) 



(31) 



/C2=0 fc™=0 



The sum over fci can be contracted by observing that fci|fci)(fci| — a'^a\ki){ki\. Then, using the commutation properties 
of the bosonic operators a, and the closure relation for the eigenstates of the harmonic oscillator, Eq. (|3l| ) takes 
the form: 



Kn^PgPT'^J dui...J du„,N{ui)...NM 

CO CO 



(32) 



k2=0 fc„,_i=0 



Repeating then the procedure shown in Eqs. (31 ),(|32|) for each index ki, we finally obtain: 

Analogously, A"^ has the form: 

^™ = rf-' Uv'i^n + l)m + ' g™(m - 1) j . 



(33) 



(34) 



Substituting now (p3[), (34) into Eqs. (ETp, (p8[) and summing over m we finally obtain Eqs. (|18|), (|19|). 
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FIG. 1. Level scheme. 
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FIG. 2. Plot of a% (solid line), and of the first term on the RHS of Eq. (19) (dashed hne) as a function of pe for Lamb-Dicke 
parameter t) = 0.1, 0.3, 0.6 and for n = Q. 
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FIG. 3. (a) Level scheme with !<;), |e), |1), stable or metastable states, |r) excited state of radiative width 7 and probability 
of decaying in the three ground states Pg, and pi, respectively. Two lasers couple |e) and |1) to |r); (b) Level scheme as in 
(a) with the addition of the excited state |2) with decay probability on |e) equal to p2, pa, respectively, P2 +P3 = 1. Two 
lasers couple |e) and |1) to |2). 
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